Coupled-channels density-matrix approach to low-energy nuclear collision dynamics: 
A technique for quantifying quantum decoherence effects on reaction observables 
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The coupled-channels density-matrix technique for nuclear reaction dynamics, which is based on 
the Liouville-von Neumann equation with Lindblad dissipative terms, is developed with the inclu- 
sion of full angular momentum couplings. It allows a quantitative study of the role and importance 
of quantum decoherence in nuclear scattering. Formulae of asymptotic observables that can reveal 
effects of quantum decoherence are given. A method for extracting energy-resolved scattering infor- 
mation from the time-dependent density matrix is introduced. As an example, model calculations 
are carried out for the low-energy collision of the 16 projectile on the 154 Sm target. 

PACS numbers: 03.65.Yz, 24.10.Eq, 24.10.-i 



I. INTRODUCTION 

Low-energy nuclear reaction dynamics has success- 
fully been treated within the stationary-state multi- 
channel scattering theory including complex potentials 
[H, 0]. However, this cannot account for quantum de- 
coherence 0-H|, which is a key aspect of irreversibility 
in open dynamical systems @, when unavoidably a lim- 
ited number of degrees of freedom and reaction chan- 
nels is used [7J. This issue can be tackled through the 
present time-dependent approach. The coupled-channels 
density-matrix (ccdm) technique was first introduced in 
studies of quantum molecular dynamics [8j], and has re- 
cently been applied to investigate the coupling-assisted 
quantum tunneling in heavy-ion fusion [H, [l(| ■ Ref . [H[ 
provides a didactic discussion on the CCDM approach. 
In contradistinction to this approach, most of the dy- 
namical models 19] of dissipative nuclear collisions do 
not treat the relative motion of the nuclei quantum- 
mechanically and/or use incoherent (statistically aver- 
aged) rather than decoherent (partially coherent) reac- 
tion channels. 

Figure Q] illustrates the innovative ideas of the CCDM 
approach. The situation can be imagined as an orches- 
tra (representing the reduced quantum system), where 
the director plays the role of the relative motion and the 
musicians correspond to a few intrinsic, low-lying collec- 
tive states. Now imagine an airplane appearing overhead, 
representing the environment of innumerable nucleonic 
states. When the airplane is not present, the orchestra 
plays a marvelous music, all musicians are in sync, like 
in a coherent quantum superposition. But when the air- 
plane approaches, the listener percieves two effects due 
to interference: (i) the music gets attenuated (dissipa- 
tion) and, most importantly, (ii) the quality of the music 
changes, as the musicians play out of sync (decoherence). 

Decoherence, which always accompanies dissipation in 
open quantum systems [6(, means dynamical dislocaliza- 
tion of coherent quantum superpositions due to entan- 
glement of the system with its environment [l2j . Co- 
herent quantum superpositions are the basis of the co- 
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FIG. 1: A low-energy nuclear collision represented by an open 
quantum system (relative motion + a few intrinsic, low-lying 
collective states). The bath of single-particle states surround- 
ing a giant resonance state represents the environment. It 
gradually destroys the coherent quantum superposition of the 
reduced-system collective states, as the nuclei approach. 



herent coupled-channels approach to near-barrier reac- 
tion dynamics, which manifest themselves through ex- 
perimental fusion-barrier distributions [llj]. Dissipation 
of energy and angular momentum is revealed in heavy- 
ion deep-inelastic scattering that occurs at near-barrier 
energies as well jlij . Cross sections of deep-inelastic col- 
lisions also indicate loss of angular-momentum coherence 
[lol Hjuf - While the coherent coupled-channels calcula- 
tions [17[ are able to explain several collision observables, 
major problems are unresolved. Foremost is the inabil- 
ity to describe the elastic and quasi-elastic scattering 
and fusion processes simultaneously [l8{ and the related, 
more recent failure to describe consistently below-barrier 
quantum tunneling and above-barrier fusion yields [l9| . 
New, precise fusion measurements have inevitably led 
to phenomenological (sometimes contradictory) adjust- 
ments [20l - [22j to stationary-state coupled-channels mod- 
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els to fit the experimental data, but without a physically 
consistent foundation. Complementary to fusion exper- 
iments, precision measurements of back-angle scattering 
energy spectra f23T - |25j and quasi-elastic barrier distribu- 
tions |26[ clearly indicate that complex, dynamical pro- 
cesses (other than the low-lying collective excitation of 
the interacting nuclei (27j ) play an important role in low- 
energy heavy- ion collisions. 

This paper suggests that quantum decoherence and 
dissipation should be simultaneously included in a consis- 
tent description of low-energy reaction dynamics, when 
a restricted set of (relevant) degrees of freedom is em- 
ployed. A possible description is the CCDM approach 
whose formalism is described in Sec. II. For the sake of 
simplicity and a specific application, without loss of gen- 
erality, the collision of an inert spherical 16 O projectile on 
a deformed 154 Sm target is studied. It includes the cou- 
plings between the relative orbital angular momentum 
L of the reactants and the spin / of the ground-state 
rotational-band levels of 154 Sm. Formulae of asymptotic 
observables that can reveal effects of quantum decoher- 
ence are obtained. A method for calculating the energy- 
resolved observables is presented. In Sec. Ill, model cal- 
culations are carried out and discussed, while a summary 
is given in Sec. IV. 

II. FORMALISM 
A. Initial density matrix 

The origin of the reference frame is in the overall cen- 
ter of mass. The vector r refers to the relative vector 
between the projectile and the target, while f and £ are 
respectively the orientation angles of r and the target 
intrinsic symmetry axis, with respect to the laboratory 
fixed z-axis defined by the beam direction. 

Initially, the target is considered to be unpolarized 
at the ground-state of spin I , being the rotational 
state the Yi m Io (0 spherical harmonic. It is cou- 
pled to the orbital motion described with Ylm l (?)- 
The total angular-momentum state, \IqL; JM) = 
T,M L M Io C lm l i m Io y lm l (r)Y IoMlQ (£), describes the an- 
gular variables, denoting C the Clebsch-Gordan coeffi- 
cients. The radial motion is described with a Gaussian 
ipk ( r ) wave-packet, centered at ro with the average ra- 
dial wave-number ko and the spatial dispersion a: 

^ ko (r) = Nexp[-^^-]e ik ^, (1) 

where the constant N is specified below. Thus, the initial 
state is \x) = J^LjM^koi^l-^oLiJM), and the initial 
density operator is po = (2/q + I) -1 |x)(x|- The latter 
can be represented as 

po = Yl \r)\a)p r J a ,(t = 0)(a'\(s\, (2) 

a, a' ,rs 



where a = (IL; JM), \a) and \r) are the coupled angular- 
momentum basis and the discrete grid-basis describing 
the internuclear separations, respectively. The initial 
density matrix is: 

p r a s a ,(t = 0) = (2/ + I)^iV 2 exp[-^?^]e lfcor 

x exp[- (5 ~ r 2 o) V ' fc ° s 5 IIo Spi , (3) 

where N is determined from the normalization condition 
V rr = I 

B. Time evolution 

1. Non-dissipative case 

The Liouville-von Neumann master equation dictates 
the time evolution of the density matrix operator p(t) 
with the initial value ([2]). It reads as 

ihf t =[H,p], (4) 

where H is the system (Hermitian) Hamiltonian speci- 
fied below. Inserting the expansion ([2]) for p(t) in Q, 
and projecting onto the basis states, we get a system of 
coupled equations for the matrix elements p T a a i(t): 

^p r ^ = E ( E % p%' <>% H K' ). ( 5 ) 

t/3 

with the initial values ([3]) . The system Hamiltonian con- 
tains different pieces: (i) the radial kinetic energy opera- 
tor T, (ii) the total bare (Coulomb+nuclear+centrifugal) 
nucleus-nucleus real potential U, (iii) the total 
(Coulomb+nuclear) real coupling potential V between 
the relative motion and the intrinsic rotational states of 
the target, and (iv) the target intrinsic Hamiltonian Hq. 
In terms of these operators, the coupled equations (JSJ) 
read as 

•rs \ y ( mrt ts rt rrits \ 

lhp aa , = 2^( T Paa< ~ Paa' T ) 

t 

+ {U a (r) - U a ,(s)]p r a s al 

+ ^2[V aP (r)p% l - p r ^V^(s)} 

P 

+ (e a - e a ,)pl s a ,. (6) 



2. Dissipative case 

In practice, however, a truncated model space of re- 
action channels (reduced system) is employed [28]. The 
impact of the excluded model space (environment) on 
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the reaction dynamics is usually treated through com- 
plex potentials 0, [H, HH , making the system Hamiltonian 
non-Hermitian. This results in loss of probability and en- 
ergy from the reduced system to the environment, but the 
quantum coherence in the basis states of the reduced sys- 
tem is preserved, as recently demonstrated in Ref. Q ■ It 
is inconsistent with the irreversible dynamics of an open 
quantum system Q, as energy dissipation always goes 
hand in hand with quantum decoherence @-|5|. 

Irreversibility can be consistently described by the 
Lindblad master equation [29H32l |. Here, quantum de- 
coherence and dissipation are incorporated into the dy- 
namics through a dissipative Liouvillian, ihCo p, which 
is added to the r.h.s. of Eq. (j4]). Cd p reads as 



for matrix elements of the reduced-system states, whereas 



E 



{c„pct--[ctc u ,p] ) 



(7) 



where [. . .] and [...]+ denote the commutator and anti- 
commutator, respectively. Each C v is a Lindblad oper- 
ator for a dissipative coupling, physically motivated ac- 
cording to the specific problem. It is assumed that each 
dissipative coupling v = {aa 1 ) between given states \a) 
and \a') has an associated (local) radial-dependent decay 
rate T^ a , Q|, i.e., the spontaneous emission Lindblad 
operator C aa ' = \/T r J a ,\a)(a'\. 

In order to also describe decay to the Hilbert space 
of excluded, intrinsic degrees of freedom, environmental 
states [34], [111 are considered in addition to the basis \a) 
of the reduced system. All these states are assumed to 
be orthonormal, resulting in an enlarged basis \a). Using 
this and the discrete grid-basis, the matrix elements of 
([7]) read as [H] 

( r p,\ rs X \ * f-prr rs Apss 

V-DPlaa' — u aa' / J y L an Pup, \/ L ay 

-^E( r ^ + r ^')^" ( 8 ) 



rr_ 

fXOt 



where runs over all the a states, and r~ = Ylu^a ^ 

37) - 

We now distinguish two kinds of density matrix ele- 
ments p™ a , , namely one type associated with the reduced- 
system states I a) only, and another type related to the 
environmental states \a), a / a. These obey the follow- 
ing equations of motion: 

A^S \ / rrirt tS rt rrttS \ 

™p aa , = 2^( T p aa > - p aa > T ) 
t 

+ [U a (r) - U a ,(s)]p™ a , 

+ ^[V a p(r) P y al -p? p V Pa ,(8)] 



p 

+ (e Q 



l P aa > 



cy / j \ flQ. 1 fJ,Ot' ) rOLOt' J J 



Paa' 



- V ( r rr - 



)PI\ 



(10) 



(9) 



for matrix elements involving the environmental states, 
i.e., either a or a' ^ a, at least. It is worth emphasizing 
that the environmental states are not reaction channels, 
but auxiliary states pjij supplying a probability drain 
and rendering the reduced-system states decoherent. 

The initial values for Eqs. (HJ) are given by ((3|), while 
for Eqs. (fT0| these are zero as the nuclei are initially 
(far apart) at the ground-states. Hence, the off-diagonal 
elements p™ a , in (fT0"|) remain zero. However, the diagonal 
terms absorb probabilities only, provided there is no 
flux back from the environment to the reduced system. 

In the calculations below, it will be considered that 
the off-diagonal elements of the decay-rate matrix L are 
nonzero only for transitions from the reduced system to 
the environment, i.e., L^, a ^ a. (Still, environment- 
induced transitions among the states of the reduced sys- 
tem may occur.) The decay rates are given by T r ^ a — 
W a (r)/h, where W a (r) > are decay functions. These 
are here treated as empirical functions, but it is hoped 
that a microsco pic theory on damping of collective ex- 
cited states |38H-4Qj . including dynamical modifications 
of excited state properties due to the close proximity of 
other nucleus [41|, can provide them. 

Various types of environments can coexist in a nuclear 
collision, which may be specific to particular degrees of 
freedom, such as isospin asymmetry or weak binding. 
Among these environments, which can be coupled to spe- 
cific states or to all states of the reduced system, are (i) 
the high level-density of one- and multi-nucleonic exci- 
tations in different mass/charge partitions (transfer), (ii) 
the continuum of non-resonant decay states of weakly- 
bound nuclei (breakup), and (iii) the innumerable nuclear 
molecular (compound nucleus) states (fusion). These 
can be treated separately, and their effects can be dis- 
tinguished within the CCDM approach. 



C. Asymptotic observables 

Having the solution of Eqs. ©, after a long period 
of time tf, when the centroid of the recoiled body of 
the wave-packet is at a large internuclear distance and 
a quasi-stationary probability current establishes in all 
the a channels, we then calculate asymptotic observables. 
For instance, these can be the angular distribution of the 
target excitations, and their integrated values. These are 
calculated as follows. 

We now introduce the projector Pim, = \IMi)(IMj\ 
associated with a specific state of the target. The new 
operator p(tf) Pim, describes the scatterred waves in this 
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target state 



Ktf)PlMr = (IMx\p(tf)\IMj 



P 



C J L Z Ml YLm(f) |r) p» x {t f ) 



x ( s l C'L'm'IMr 



(11) 



where p = (r, s, L, m, J, M, L', m', J', M'), 7 = (IL; JM) 
and A = {IL'\J'M'), The reduced density matrix, 
pZ\(tf), is normalized with its trace, i.e., *^2 ri p^(tf). 

We note that is still an operator in r and s. The 
radial projector P r > = \r') (r'\ associated with a specific 
separation between the nuclei is then introduced. With 
this projector we now act on (fTTj) and set f — s — r' . A 
new operator is thus obtained, whose partial trace (sum 
over all separations r') provides the probability for pro- 
ducing the target in state (I, Mi) with the relative coor- 
dinate in the direction f'\ 
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X C-L'm'IM! Y£> m >(f), 



(12) 



where q = (L,m,J,M,L',m',J\M') and 5 7 a(*/) = 
Sr' P7A The latter contains information about the 

coherence of angular momenta. 

Integrating (|T2l over all directions f' of solid angles, 
and summing over all Mi, the total probability for pro- 
ducing the target in state I (population) is obtained: 



(13) 



Mj LrnJM 



D. Energy-resolved observables 

The observables (fT2")l and (|13[) correspond to average 
values for the range of energies contained in the incident 
wave-packet JT]). The energy-resolved scattering infor- 
mation can be obtained using a window operator (42j. 
The key idea is to calculate, for definite 7A indices, the 
energy spectrum V{Ek) of the initial and final reduced 
density matrices. Ek is the centroid of a total energy bin 
of width 2e. A matrix of reflection coefficients, TZ^\(Ek), 
is determined by the ratio 



n lX {E k ) 



v. 



final 
7 A 



E 7 V™ Ual (E k ) 



(14) 



which replaces the matrix S 7 \(tf) in (fT2|) and (p~3]l . Ex- 
pression (|14p generalizes the wave-packet formulation of 
the reflection coefficient [43|, HH . 

The energy spectrum V(Ek) = Tr(Ap), where Tr de- 
notes a partial trace involving the radial indices only, and 
A is the window operator j42| : 



A(E k ,n,e) = 



(15) 



where H is the system asymptotic Hamiltonian, and n 
determines the shape of the window function. As n is 
increased, this shape rapidly becomes rectangular with 
very little overlap between adjacent energy bins (42|, re- 
maining the bin width constant at 2e. The spectrum is 
constructed for a set of Ek where Ek+i — E k + 2e. Thus, 
scattering information over a range of incident energies 
can be extracted from a time-dependent density matrix 
that has been calculated on a grid. 



1. Example 

Figure [5] shows for a single-channel CCDM calculation 
with L = @: (a) the energy spectrum (n = 4, 2e = 1 
MeV) of the initial (solid line) and final (dashed line) 
density matrices, and (b) the final-to-initial ratio of the 
energy spectrum providing energy-resolved reflection co- 
efficients (full squares) . These very well agree with those 
(dotted line) of a time-independent Schrodinger equation 
with a short-range imaginary potential or an ingoing- 
wave boundary condition at small radii, as implemented 
in the CCFULL code |45|. 

In the CCDM calculation above, the grid (r = — 100 
fm) was evenly spaced with 512 radial points. The in- 
cident wave-packet was initially centered at ro — 50 fm, 
with width a = 5 fm, and was boosted toward the poten- 



tial with the appropriate average kinetic energy for the 
total energy E = 60 MeV. 

The time evolution of the density matrix was carried 
out using a Faber propagator [46| and the Fourier method 
[47| for the commutator between the kinetic energy and 
density operator. The time step for the density-matrix 
propagation was At = 10~ 22 s. 

The form of the bare nuclear potential between the 
16 and 154 Sm nuclei is a Woods-Saxon potential with 
(V ,r ,a ) = (-165 MeV, 0.95 fm, 1.05 fm). The 
Coulomb potential is that for two point charges. These 
yield a s-wave Coulomb barrier of Vb = 59.41 MeV at 
the radius Rb = 10.81 fm. 

The irreversible capture of the two nuclei in the 
nucleus- nucleus potential pocket inside the Coulomb bar- 
rier(fusion) is caused by an environmental coupling 
[9lll0l| with a decay function W(r) taken as a Fermi func- 
tion with depth 10 MeV and diffuseness 0.1 fm, located 
at the pocket radius of 5.77 fm. 



III. MODEL CALCULATIONS 

In the model calculations the 16 O projectile was taken 
to be inert and the 154 Sm target was allowed to be excited 
up to the 4 + state of the ground-state rotational band. 
The all order nuclear coupling of the ground-state + to 
the states 2 + and 4 + , with excitation energies E 2 + = 0.08 
MeV and E^ = 0.27 MeV, has a macroscopic deformed 
Woods-Saxon form with the radius parameter of 1.06 
fm, and the deformation parameters of — 0.322 and 
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FIG. 2: (a) Energy spectrum of the initial (solid line) and 
final (dashed line) density matrices for the elastic scattering 
of 16 off the 154 Sm target at the average total energy Eq = 60 
MeV. (b) Energy-resolved reflection coefficients provided by 
the final-to-initial ratio of the energy spectrum (full squares). 
These very well agree with those of a stationary Schrodinger 
equation (dotted line). 



@4 = 0.027. The Coulomb coupling includes terms up to 
second order with respect to and to the first order of 
/?4. The total coupling-potential matrix in the coupled 
angular-momentum basis is calculated using the CCFULL 
code [45j and the Kermode- Rowley matrix technique [48j . 

With such a coupling Hamiltonian, the time propaga- 
tion on the grid employed in Fig. [5] is very time consum- 
ing and memory demanding, beyond the present limit of 
our computational capability. That is why the present 
calculations aim at exploring qualitative effects only, for 
which a smaller grid (r = — 40 fm) with 64 evenly 
spaced radial points suffices. The initial wave-packet is 
then centered at ro = 25 fm, with width a = 3 fm and 
the incident, average total energy E = 60 MeV. Relative 
partial waves up to 20H are included in the calculations. 

Two types of calculations (without energy projection) 
are carried out including: (i) only the effects of the fu- 
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FIG. 3: Time evolution of the measure of coherence in the re- 
duced system. While the fusion environment essentially pre- 
serves coherence (solid line), it is destroyed by the surface 
environment (dashed line). 



sion environment highlighted above, and (ii) in addi- 
tion the effects of a surface environment specified below. 
While the fusion environment is coupled to all the 154 Sm 
states, the surface environment is considered to be cou- 
pled to the ground-state only. The latter can be associ- 
ated with complex, multi-nucleon/cluster transfers from 
the ground-state of the colliding nuclei to other mass (or 
charge) partitions. The corresponding decay function is 
taken as a Gaussian with width of 1 fm, centered at the 
contact radius that is estimated as 1.2(I6 1 / 3 + 154 1 / 3 ) 
fm. This function is physically motivated by the spatial 
localization of transfer processes in heavy-ion reactions 
[49l j . The measure of coherence 0, H(| in the reduced 
system is the ratio Tr(p 2 )/[Tr(p)] 2 , whose time evolution 
is presented in Fig. [3] The fusion environment essen- 
tially preserves coherence (solid line), while the surface 
environment results in decoherence (dashed line). Com- 
paring the calculation (ii) to (i), we learn how the surface 
environment-induced decoherence impacts on the asymp- 
totic observables (fT2")l and (|I3I) . 

Figure @] shows the angular distribution of the 
154 Sm excitations, which corresponds to expression (fT2|) 
summed over all Mj. The solid and dashed lines are re- 
sults of the calculations (i) and (ii) , respectively. Clearly, 
the surface environment varies the quantum interference 
effects, destroying the coherence of relative partial waves 
and changing by a few degrees the minimum of the in- 
elastic probability distributions [Figs. 0Jb) and (c)]. It 
also affects significantly the asymptotic population of the 
154 Sm states and the fusion probability, as presented in 
Table |TJ The surface environment-induced decoherence 
hinders the probability flow from the elastic to the in- 
elastic and fusion channels. 
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FIG. 4: Angular distribution of the 154 Sm excitation probabil- 
ities, when an inert ie O projectile scatters off the 154 Sm target 
at the average total energy Eo = 60 MeV: (a) + , (b) 2 + , and 
(c) 4 + . The solid lines are outcomes including the effects of 
the fusion environment only, while the dashed lines include 
in addition the effects of the surface environment. The latter 
destroys the coherence of relative partial waves, changing by a 
few degrees the minimum of the inelastic distributions (small 
figures inserted). 



TABLE I: Asymptotic population |I3| of the 154 Sm states 
and the fusion probability for calculations including (i) only 
the fusion environment (Environment I), and (ii) in addition 
the surface environment (Environment 2) . The latter hinders 
the probability flow from the elastic to the inelastic and fusion 
channels, due to decoherence. 





Environment I 


Environments 1 and 2 


States 


Population 


Population 


0+ 


0.8376 


0.8770 


2+ 


0.1206 


0.0961 


4+ 


0.0418 


0.0268 


Fusion Probability 


3.127E-3 


9.106E-4 



IV. SUMMARY 

The innovative, coupled-channels density-matrix ap- 
proach to low-energy reaction dynamics has been de- 
veloped further, including full angular momentum cou- 
plings and a method for extracting energy-resolved scat- 
tering information from the time-dependent density ma- 
trix. These developments have enabled a first exploration 
of decoherence effects on asymptotic observables, such as 
the angular distribution of the target excitations, and 
their integrated values. These are significantly affected 
by decoherence induced by a surface environment (re- 
lated to complex, multi-nucleon/cluster transfers), which 
changes by a few degrees the minimum of the back- angle 
inelastic probability distributions and hinders the prob- 
ability flow from the elastic to the inelastic and fusion 
channels. To test the present theory against experiments, 
both extensive calculations (with energy projection) and 
high-precision measurements of fusion and scattering ob- 
servables (including those investigated here) are required. 
It is hoped that decoherence effects can help resolve ma- 
jor problems in low-energy nuclear reaction physics, such 
as the current inability to simultaneously describe fusion 
and scattering measurements of heavy ions within the 
standard, coupled-channels framework. 



Acknowledgments 

The author thanks Jeff Tostevin for discussions, and 
Ron Johnson for constructive comments on the impor- 
tance of the density-matrix off-diagonal elements (with 
respect to the total angular momentum) for calculating 
the angular distributions. Useful discussions with partic- 
ipants in the ECT* Workshop on Decoherence in Quan- 
tum Dynamical Systems, Trento, April 26-30, 2010, are 
also acknowledged. The work was supported by the UK 
Science and Technology Facilities Council (STFC) Grant 
No. ST/F012012/1. 



7 



[1] H. Feshbach, Ann. Phys. (N.Y.) 19, 287 (1962), and ref- 
erences therein. 

[2] J.R. Taylor, Scattering Theory (John Wiley and Sons, 
Inc., 1972). 

[3] E. Joos, H.D. Zeh, C. Kiefer, D. Giulini, J. Kupsch and 
1.0. Stamatescu, Decoherence and the Appearance of a 
Classical World in Quantum Theory (Springer, Heidel- 
berg, 2003). 

[4] W.H. Zurek, Rev. Mod. Phys. 75, 715 (2003). 

[5] M. Schlosshauer, Decoherence and the Quantum-to- 

Classical Transition (Springer, Berlin, 2007). 
[6] H.-P. Breuer and F. Petruccione, The Theory of Open 

Quantum Systems (Oxford University Press, New York, 

2002). 

[7] A. Diaz- Torres, Phys. Rev. C 81, 041603(R) (2010). 
[8] L. Pesce and P. Saalfrank, Chem. Phys. 219, 43 (1997); 

J. Chem. Phys. 108, 3045 (1998). 
[9] A. Diaz- Torres et al., Phys. Rev. C 78, 064604 (2008). 
[10] A. Diaz- Torres et al., AIP Conf. Proc. 1098, 44 (2009). 
[11] A. Diaz- Torres, larXiv:1009.0520l [nucl-th], 
[12] H.D. Zeh, Seminaire Poincare 1, 115 (2005). 
[13] M. Dasgupta et al., Annu. Rev. Nucl. Part. Sci. 48, 401 
(1998). 

[14] D.J. Hinde et al., Nucl. Phys. A 834, 117c (2010), and 

references therein. 
[15] A.Y. Abul-Magd and M.H. Simbel, Phys. Lett. B 83, 27 

(1979). 

[16] A. De Rosa et al., Phys. Rev. C 40, 627 (1989). 
[17] A.B. Balantekin and N. Takigawa, Rev. Mod. Phys. 70, 
77 (1998). 

[18] J.O. Newton et al., Phys. Rev. C 70, 024605 (2004). 
[19] M. Dasgupta et al., Phys. Rev. Lett. 99, 192701 (2007), 

and references therein. 
[20] I. Boztosun and W.D.M. Rae, Phys. Lett. B 518, 229 

(2001). 

[21] S. Misicu and H. Esbensen, Phys. Rev. Lett. 96, 112701 
(2006). 

[22] T. Ichikawa, K. Hagino and A. Iwamoto, Phys. Rev. Lett. 

103, 202701 (2009). 
[23] L.R. Casques et al., Phys. Rev. C 76, 024612 (2007). 
[24] M. Evers et al., Phys. Rev. C 78, 034614 (2008). 



[25] C.J. Lin et al., Phys. Rev. C 79, 064603 (2009). 
[26] E. Piasecki et al., Phys. Rev. C 80, 054613 (2009). 
[27] S. Yusa, K. Hagino and N. Rowley, Phys. Rev. C 82, 
024606 (2010). 

[28] S. Karataglidis and K. Amos, Phys. Lett. B 660, 428 
(2008). 

[29] G. Lindblad, Comm. Math. Phys. 48, 119 (1976). 

[30] V. Gorini et al., J. Math. Phys. 17, 821 (1976). 

[31] A. Sandulescu, H. Scutaru and W. Scheid, J. Phys. A 20, 

2121 (1987); Ann. Phys. (N.Y.) 173, 277 (1987). 
[32] R. Kosloff, M.A. Ratner and W.B. Davis, J. Chem. Phys. 

106, 7036 (1997). 
[33] C. Scheurer and P. Saalfrank, J. Chem. Phys. 104, 2869 

(1996). 

[34] I. Burghardt, J. Phys. Chem. A 102, 4192 (1998). 

[35] R.A. Bertlmann et al., Phys. Rev. A 73, 054101 (2006). 

[36] Expression (6) in Ref. |{j] is incorrect, but the correction 

of this error does not affect any conclusions of that paper. 
[37] C. Scheurer and P. Saalfrank, Chem. Phys. Lett. 245, 

201 (1995). 

[38] M. Matsuo et al., Nucl. Phys. A 649, 379c (1999). 
[39] L. Guo et al., Nucl. Phys. A 753, 136 (2005). 
[40] E. Litvinova et al., Phys. Rev. C 75, 064308 (2007). 
[41] A.B. Mcintosh et al., Phys. Rev. Lett. 99, 132701 (2007). 
[42] K.J. Schafer and K.C. Kulander, Phys. Rev. A 42, 5794 

(1990); Comp. Phys. Comm. 63, 427 (1991). 
[43] D.J. Tannor, Introduction to Quantum Mechanics: A 

Time-Dependent Perspective (University Science Books, 

Saulito, 2007). 
[44] K. Yabana, Prog. Theor. Phys. 97, 437 (1997). 
[45] K. Hagino, N. Rowley and A.T. Kruppa, Comp. Phys. 

Comm. 123, 143 (1999). 
[46] W. Huisinga et al., J. Chem. Phys. 110, 5538 (1999). 
[47] R. Kosloff, Annu. Rev. Phys. Chem. 45, 145 (1994). 
[48] M.W. Kermode and N. Rowley, Phys. Rev. C 48, 2326 

(1993). 

[49] G.R. Satchler, Direct Nuclear Reactions (Oxford Univer- 
sity Press, Oxford, 1983) p. 684. 

[50] K. Blum, Density Matrix Theory and Applications (Sec- 
ond Edition, Plenum Press, New York, 1996) p. 39. 



